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ABSTRACT 

Direct A-body simulations of star clusters in a realistic Milky Way-like potential are carried out 
using the code NBDDY6. Based on these simulations a new relationship between scale size and galac- 
tocentric distance is derived: the scale size of star clusters is proportional to the hyperbolic tangent 
of the galactocentric distance. The half-mass radius of star clusters increases systematically with 
galactocentric distance but levels off when star clusters orbit the galaxy beyond ~40 kpc. These sim- 
ulations show that the half-mass radius of individual star clusters varies significantly as they evolve 
over a Hubble time, more so for clusters with shorter relaxation times, and remains constant through 
several relaxation times only in certain situations when expansion driven by the internal dynamics of 
the star cluster and the influence of the host galaxy tidal field balance each other. Indeed, the radius of 
a star cluster evolving within the inner 20 kpc of a realistic galactic gravitational potential is severely 
truncated by tidal interactions and does not remain constant over a Hubble time. Furthermore, the 
half-mass radius of star clusters measured with present day observations bears no memory of the 
original cluster size. Stellar evolution and tidal stripping are the two competing physical mechanisms 
that determine the present day size of globular clusters. These simulations also show that extended 
star clusters can form at large galactocentric distances while remaining fully bound to the host galaxy. 
There is thus no need to invoke accretion from an external galaxy to explain the presence of extended 
clusters at large galactocentric distances in a Milky Way-type galaxy. 

Subject headings: Galaxy: star clusters: general: galaxies: star clusters - galaxies: dwarf - galaxies- 
stars: evolution 



1. INTRODUCTION 

Star clusters are an increasingly diverse family. During 
the last decade the discovery of stellar systems brighter, 
larger, and more massive than the "standard" star clus- 
ter has blurred the distinction between globular clusters 
and dwarf galaxies. Ultra-compact dwarfs (UCDs) are 
the prime example of stellar systems with physical pa- 
rameters between those of globular clusters and dwarf 
elliptical galaxies (Hilker et al. 1999; Drinkwater et al. 
2000). At the other end of the luminosity range, faint and 
extended star clusters have also been discovered, notably 
as satellites of the Andromeda galaxy (Huxor et al. 2005; 
Mackey et al. 2006) and NGC 1023 (Larsen & Brodie 
2000). 

These newly discovered stellar systems have bridged 
the gap in physical size thought to exist between globu- 
lar clusters and compact elliptical galaxies (Gilmore et al. 
2007). The well defined linear relations between physical 
size and total magnitude for galaxies with masses greater 
than 10 8 Mo can be extrapolated to UCDs. However, 
dwarf galaxies and star clusters form two branches in 
the size-magnitude plane where these two physical pa- 
rameters are uncorrelated (Misgeld & Hilker 2011 their 
Figure 1; see also McLaughlin 2000). 

In recent years, studies based on the relatively wide 
field of view and superb resolution of the Advanced 
Camera for Surveys onboard the Hubble Space Telescope 
(HST), have found a strikingly constant median effective 
radius, or cquivalently half-light radius, for extragalac- 
tic globular clusters of < >~ 3 pc. Jordan et al. 
(2005) provide a prime example of such work as they ac- 
curately determined the structural parameters of thou- 
sands of globular clusters associated with 100 early type 



galaxies of the Virgo cluster (ACS Virgo Cluster Survey). 
Masters et al. (2010) replicated this work for 43 galaxies 
in the Fornax cluster (ACS Fornax Cluster Survey). One 
of the main findings of these papers is that thousands of 
star clusters spanning more than four magnitudes in lu- 
minosity have the same median value of < >^ 3 pc. 

Why is the size distribution of star clusters narrowly 
centered around three parsecs? Why do only a few clus- 
ters become extended with effective radii of ten parsecs 
or more? In this work, advanced A-body models are car- 
ried out with the aim of determining the most important 
physical mechanisms that mold the characteristic radii 
of star clusters. The impact of the host galaxy tidal field 
on the size of orbiting star clusters is probed in detail by 
evolving several models at different galactocentric dis- 
tances (Rqc)- 

The empirical qualitative dependence between size and 
galactocentric distance of star clusters has been clearly 
established in several observational studies beginning 
with the work of Hodge (1960, 1962). The A-body mod- 
els that have been performed allow us to quantify the 
influence of the tidal field, generated by a Milky- Way 
or M31 type galaxy, on satellite star clusters. We thus 
determine a new relation between the scale sizes of star 
clusters and galactocentric distance. This relation is a 
proxy for the host galaxy gravitational potential. 

2. THE MODELS 

All A-body simulations are performed using the 
NB0DY6 code (Aarseth 1999; Aarseth 2003). This code 
performs a direct integration of the equations of motion 
for all N stars and binaries in a star cluster and includes 
a comprehensive treatment of stellar evolution (Hurley et 
al. 2000, 2005). This code also includes a detailed han- 
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dling of binaries accounting for close encounters, mergers, 
and the formation of three and four body systems (Tout 
2008; Hurley 2008; Mardling 2008; Mikkola 2008). 

The version of NB0DY6 that is used in this work was 
specially modified to run on a Graphic Processing Unit 
(GPU: Nitadori k Aarseth 2012). During the last decade 
the clock rates of Central Processing Units (CPUs) have 
been practically stagnant while CPUs provide a proven 
alternative for high-performance computing (Barsdell et 
al. 2010), and particularly for iV-body codes (Hamada 
et al. 2009). The models are run at the Center for As- 
trophysics and Supercomputing of Swinburne University. 
The GPUs in use are NVIDIA Tesla S1070 cards. Earlier 
versions of this code (i.e. NB0DY4) were run on special- 
purpose GRAPE hardware (Makino et al. 2003) but the 
performance of the GPU version of this code is compara- 
ble or superior to previous efforts to improve computing 
time. The calculations are carried out in ./V-body units, 
i.e. G = 1 and — 4Eo = 1, where Eq is the initial energy 
(Heggie & Mathieu 1986). The results are scaled back 
to physical units once the computation of the models is 
complete. 

2.1. Numerical Simulations Set Up 

The initial set up for the simulations carried out here is 
similar to the work of Hurley & Mackey (2010). All simu- 
lations have an initial number of particles of N — 10 5 . Of 
these, 5% are primordial binary systems, that is, 95 000 
single stars and 5000 binary systems. The most massive 
star has a mass of M max = 50M Q while the minimum 
mass for a star is M min = 0.1M Q . The initial mass 
distribution for all stars follows the stellar Initial Mass 
Function (IMF) of Kroupa et al. (1993). The models 
start off with a total initial mass of M to t ~ 6.3 x 10 4 M Q . 
The N = 10 5 stellar systems have the initial spatial dis- 
tribution of a Plummer sphere (Plummer 1911) and an 
initial velocity distribution that assumes virial equilib- 
rium. The modeled clusters stars have a metallicity of 
Z = 0.001 or [Fe/H] « -1.3. 

These simulations start after the clusters will have un- 
dergone expansion as the result of the removal of resid- 
ual gas left over from star formation. Thus, at t = all 
stars are assumed to have formed and be on the zero age 
main sequence, with no gas present. The initial three- 
dimensional half-mass radius is 6.2 pc. 

The initial set up for each simulation is exactly the 
same with the exception of the initial galactocentric dis- 
tance. Individual simulations of star clusters evolving on 
circular orbits at different galactocentric distances were 
obtained, i.e. R G c = 4, 6, 8, 8.5, 10, 20, 50 and 100 
kpc. The initial plane of motion of the star clusters is 
22.5 degrees from the plane defined by the disk. 

2.2. Galactic Tidal Field Model 

Earlier models from A r -body simulations have made 
the simplifying assumption that the gravitational poten- 
tial of the host galaxy can be represented by a central 
point mass only (e.g. Vesperini & Heggie 1997; Baum- 
gardt 2001; Hurley & Bekki 2008). The version of 
NBDDY6 used here models the bulge as a point source 
but also includes a halo and a disk as components of the 
Milky- Way galaxy and its gravitational potential. This 
more realistic implementation of the host galaxy poten- 



tial has been incorporated in recent work (e.g. Kiipper 
et al. 2011). 

To model the disk NBDDY6 follows Miyamoto & Na- 
gai (1975) whom combined the potential of a spherical 
system (Plummer 1911) and the potential of a disk-like 
mass distribution (Toomre 1963) into a generalized ana- 
lytical function that elegantly describes the disk of spiral 
galaxies: 



$(r, z) = 



GM 



v /r 2 + [a+ v /"(z 2 + 6 2 )] 2 



(1) 



Here a is the disk scale length, b is the disk scale height, 
G is the gravitational constant, and M is the mass of the 
disk component. The elegance of this formulation resides 
in the fact that it can model both an infinitely thin disk 
or a sphere by varying the two scale factors a and b. The 
values used here are a = 4 kpc for the disk scale length, 
and b = 0.5 kpc for the vertical length (Read et al. 2006). 
Formally, the Miyamoto & Nagai disk extends to infinity. 
However, the strength of the disk potential asymptotes 
towards zero at large radii: with a = 4 kpc and b = 0.5 
kpc, the density at 40 kpc drops to 0.1% of the central 
value. 

The masses of the bulge and disk are 1.5 x 1O 1O M0 and 
5 x 10 10 M Q , respectively (Xue et al. 2008). The galac- 
tic halo is modeled by a logarithmic potential. When 
combined together, the potential of the halo, disk, and 
bulge are constrained to give a rotational velocity of 220 
km/s at 8.5 kpc from the galactic center (Aarseth 2003). 
A detailed discussion of the galactic tidal field model is 
given by Praagman et al. (2010). 

A tidal radius for a star cluster in a circular orbit about 
a point-mass galaxy, initially postulated by von Hoerner 
(1957), can be approximated by the King (1962) formu- 
lation: 



( M c \ 1/3 R 



(2) 



where Mc is the mass of the cluster, and Mq is the 
mass of the galaxy. The equivalent expression for a cir- 
cular orbit in the Milky- Way potential described above 



is: 



r t ^ 



(GM C \ 1/3 
\ 2ft 2 J 



(3) 



where ft is the angular velocity of the cluster (Kiipper 
et al. 2010a). The tidal radius is where a star will feel an 
equal gravitational pull towards the cluster and towards 
the galaxy center in the opposite direction. A detailed 
study of the tidal radius of a star cluster for different 
galaxy potentials is given by Renaud et al. (2011). 

In the ./V-body code the gravitational forces owing to 
both the cluster stars and the galaxy potential are taken 
into account for all stars in the simulation, independently 
of the definition of the tidal radius. However it is com- 
mon practice to define an escape radius, beyond which 
stars are deemed to be no longer significant in terms of 
the cluster potential and thus have no further input on 
the cluster evolution. Stars are removed from the simu- 
lation when they fulfill two conditions: (i) their distance 
from the cluster center exceeds the escape radius, and (ii) 
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they have positive energy when the external field contri- 
bution is taken into account. 

In our simulations we have estimated the escape radius 
as twice the tidal radius given by Eq. (2) with Mq — 
10 11 Mq. At all times this gives a value in excess of the 
tidal radius given by Eq. (3). When presenting results, 
we only consider stars that lie within the tidal radius 
given by Eq. (3). 

3. PREVIOUS WORK 

A pioneering study of this subject was carried out by 
Vesperini & Heggie (1997) who investigated the effects 
of dynamical evolution on the mass function of globular 
clusters through simulations based on NBDDY4. This was 
the most up-to-date version of the code at the time of 
their work and only included a basic treatment of stellar 
evolution. Also, the total number of particles that Ves- 
perini & Heggie (1997) were able to simulate was limited 
to N — 4096 due to restrictions imposed by the hard- 
ware. Under the conditions used here a simulated cluster 
with only N — 4096 orbiting at Rgc— 8 kpc is dissolved 
in 1.1 Gyr. Despite the computational limitations to 
which Vesperini & Heggie (1997) were subject, they set 
an inspiring precedent for this work. Particularly rele- 
vant is the trend of increasing mass-loss with decreasing 
galactocentric distance, albeit with the galaxy potential 
modeled as a point-mass. 

Baumgardt & Makino (2003) studied the stellar mass 
function of star clusters using NB0DY4 but with more re- 
alistic particle numbers, going up to N = 131072. They 
determined that stellar evolution accounts for ~ 1/3 of 
the total mass lost by the cluster. The upper mass limit 
of the stellar initial mass function used by Baumgardt 
& Makino (2003) was M max — 15M Q , while in the sim- 
ulations presented here M max = 50M Q . Baumgardt & 
Makino (2003) found that owing to mass segregation low- 
mass stars are prone to be depleted from the star clus- 
ter. This depletion of low mass stars inverts the slope of 
the IMF. Unlike Baumgardt & Makino (2003) the mod- 
els presented here include a population of 5000 primor- 
dial binary systems. For the original spatial distribution 
Baumgardt & Makino (2003) used a King profile while 
here a Plummer sphere is used. 

4. TOTAL MASS 

NB0DY6 readily yields fundamental numerical param- 
eters of the simulated star clusters such as the number 
of stars, number of binaries, total mass, core mass, half- 
mass radius, relaxation time, and velocity dispersion, to 
name a few. 

Star clusters evolving in a galactic potential lose mass 
due to stellar evolution, two body relaxation, few-body 
encounters, and tidal interactions with the host galaxy. 
In the following sections the two main phases of mass-loss 
are discussed. 

4.1. Initial Stellar Evolution 

During the first Gyr of evolution star clusters lose large 
amounts of mass owing to stellar evolution. Massive stars 
with M~ 10M Q and above such as OB stars, luminous 
blue variables, and Wolf-Rayet stars have large mass- 
loss rates of the order of M <~ 10~ 5 M Q per year (Van- 
beveren et al. 1998; Shara et al. 2009). The mass loss 



rates of these stars is the dominant factor of their evo- 
lution. In fact, mass-loss determines whether a massive 
star becomes a supernovae or a long duration gamma-ray 
bust (Vink et al. 2000, Shara et al. 2009). Massive stars 
are short lived, for instance Wolf-Rayet stars with initial 
masses greater than 2OM are expected to survive only 
10 Myr (Hurley et al. 2000). 

The iV-body simulations plotted in Figure 1 show that 
star clusters, independently of galactocentric distance, 
lose <~ 1/3 of their mass during the first Gyr of evolution. 
This is in agreement with the findings of Baumgardt & 
Makino (2003). The stellar evolution of massive stars dis- 
cussed above is responsible for the mass loss that triggers 
an expansion of the cluster due to a reduced gravitational 
pull. The impact of mass loss on the scale size of a star 
cluster is discussed in the next section. 

4.2. Two phases of mass-loss 

The two main regimes of mass loss evident in these 
simulations are mass loss from stellar evolution and mass 
loss due to tidal stripping. The former initially increases 
the cluster size, while the latter decreases the cluster size. 
Whether or not either of the two mechanisms dominates 
will determine the size of the star cluster (Gieles et al. 
2011). 

As mentioned above, stellar evolution mass-loss occurs 
on a rapid timescale at early times. It subsequently slows 
owing to the longer evolution timescales of low-mass stars 
but maintains a steady presence throughout the ongoing 
life of a cluster. Mass-loss resulting from tidal stripping is 
linked to the two-body relaxation timescale of the cluster: 
two-body interactions gradually increase the velocities of 
the low-mass stars, pushing them to the outer regions 
where they are lost across the tidal boundary in a pro- 
cess that is often called evaporation (McLaughlin & Fall 
2008). For star clusters evolving at small galactocentric 
distances (e.g. simulations at 4 and 6 kpc) the dominant 
mechanism of mass loss is tidal stripping. These clusters 
have relatively short two-body relaxation timescales. In 
contrast, for clusters evolving at large galactocentric dis- 
tances (with longer relaxation timescales) the dominant 
mass loss mechanism is stellar evolution. It is clear from 
Figure 1 that the closer a star cluster orbits to the cen- 
ter of the galaxy, the more accentuated its mass loss is. 
Figure 1 shows how the mass-loss rate evolves towards 
an asymptotic linear behavior after a Hubble time, in 
agreement with Baumgardt & Makino (2003). Few-body 
interactions can also lead to the loss of stars through ejec- 
tions. This process can cause small-scale fluctuations in 
cluster size on short timescales. 

5. CHARACTERISTIC RADII 

Different characteristic radii are commonly used for 
star clusters. Observers obtain two dimensional half- 
light radii to high accuracy, but are limited in their choice 
of characteristic radii. For simulated clusters, different 
scale size values are readily available such as the half- 
mass radius, r/j TO , or core radius. Radii of simulated 
clusters can be expressed in two or three dimensions, the 
latter is adopted in this work. Hurley & Mackey (2010) 
show that the value of the half- mass radius derived with 
NB0DY6 is on average 1.6 times the value of an observed 
half-light radius. Moreover, different characteristic radii 
describing the core of star clusters are known to evolve 
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Fig. 1. — Mass loss of simulated star clusters at different galac- 
tocentric distances. Note how the simulated star cluster evolving 
at 4 kpc has lost all of its mass after 8.4 Gyr. 

in a self-similar fashion (Baumgardt et al. 2004). 

The definition of core radius is different in observa- 
tional and theoretical studies. In observational studies 
the core radius is generally defined as the radius where 
the surface brightness falls to half its central value (King 
1962). In our simulations, the core radius is a density- 
weighted average distance of each star to the point of 
highest stellar density within the cluster (Aarseth 2003). 

5.1. Observed and primordial half-mass radius 

The evolution of the half-mass radius of simulated star 
clusters with time is plotted in Figure 2. Clusters start 
with an initial half- mass radius of ~ 6.2 pc and undergo 
an expansion triggered by stellar evolution within the 
first Gyr as discussed above. An initial half-mass radius 
of ~ 6.2 pc might seem larger than average but after a 
Hubble time of evolution King (1962) models fitted to 
the light profiles of our models yield two dimensional 
half-light radii that are consistent with observations, i.e. 
around 3 pc (Sippel et al. 2012). 

In their original work, Spitzer & Thuan (1972) found 
that the effective radius of an isolated star cluster re- 
mains constant through several relaxation times. Their 
important result is confirmed in Figure 2 but only for the 
model evolving at 20 kpc from the galactic center. In this 
particular model the star cluster experiences only weak 
tidal interactions with the host galaxy. After the ini- 
tial phase of rapid expansion driven by stellar evolution, 
the longer-term steady expansion driven by the internal 
dynamics is balanced by the presence of the tidal field. 
Thus, the half-mass radius of this simulated star cluster 
remains constant over the last eight Gyr of evolution in 
agreement with Spitzer & Thuan (1972). 

More importantly, the effective radii of all star clusters 
evolving at Rgc < 20 kpc are significantly affected by 
tidal interactions that truncate the cluster size during 
its evolution. The results of Spitzer & Thuan (1972) are 
often misquoted in observational studies as proof that the 
effective radius remains constant over many relaxation 
times and that the Th measured today is a faithful tracer 
of the original size of the proto-cluster cloud. 



Fig. 2. — Evolution of the 3D half-mass radius of simulated 
star clusters at different galactocentric distances. A simulated star 
cluster evolving at 4 kpc from the galactic center dissolves before 
a Hubble time. Simulations of star clusters orbiting the galaxy at 
less than 20 kpc experience a truncation of their half-mass radius 
due to tidal interactions. 

As shown in Figure 2, the half- mass radii of simulated 
clusters evolving at 50 and 100 kpc are always increasing. 
For these two simulations, in virtual tidal isolation and 
with long relaxation times, their half-mass radii reach 
about twice the value of the initial half-mass radius, i.e. 
I'hm/i'hmO ~ 2 after a Hubble time (Gieles et al. 2010). 

Most observational studies of extragalactic star clus- 
ters aimed at determining their scale sizes are carried 
out with the Hubble Space Telescope due to its unique 
resolution (e.g. the ACS Virgo Cluster Survey by Jordan 
et al. 2005). The drawback of space-based detectors is 
the small field of view that only covers, in general, a 
physical scale of a few kiloparsecs in radius surrounding 
the core of the host galaxy. As shown in Figure 2 the 
inner ten kiloparsercs of a galaxy is where the tidal field 
has the strongest impact on the size of the satellite star 
cluster, with a dependence on galaxy mass. The effective 
radius of a star cluster evolving in a realistic galaxy po- 
tential loses any memory of its earlier values. As shown 
in Figure 2, only the simulated star cluster evolving at 20 
kpc has an effective radius that remains constant through 
several relaxation times which is a result of its particu- 
lar circumstances, being in a position within the galaxy 
where internally driven expansion is balanced by the ex- 
ternal truncation of the tidal field. The exact position 
at which this occurs will be dependent on the strength 
of the tidal field, and thus the galaxy model, and on the 
initial mass of the star cluster. 

In the simulations presented here a star cluster orbit- 
ing the Galactic center at 50 kpc or more is free from 
disk shocking, given that the disk has an extent of ^40 
kpc. Disk shocking has demonstrated wounding effects 
on cluster stability (Gncdin & Ostriker 1997; Vesperini 
& Hcggic 1997). In an effort to quantify the effect of 
disk shocking a full simulation was carried out with a 
star cluster in an orbit at 10 kpc from the galactic cen- 
ter where the mass of the disk was placed in a central 
spheroid component instead of a disk. This showed that 
the presence of a disk will enhance the mass loss of the 
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simulated cluster by 2.4x 10 3 M Q (13% of the total mass) 
over a Hubble time, and will make its half-mass radius 
smaller by 0.6 pc (10% of its size). Thus at 10 kpc from 
the galactic center the effect of the disk in a star cluster, 
while subtle, is clearly measurable. 

From the output of NB0DY6 we compute the half-mass 
relaxation time as: 

0.147V /TjT" fA . 

where N is the number of stars, n im the half-mass 
radius, and A = 0.4iV the argument of the Coulomb 
logarithm (Spitzer & Hart 1971; Binney & Tremainc 
1987). Note that Giersz & Heggie (1994) found a value 
of A = 0.11 AT. The evolution of the relaxation time for 
each star cluster is homologous to the evolution of the 
half-mass radius with time shown in Figure 2. For all 
clusters it reaches a maximum value at the point when 
expansion driven by mass loss is equivalent to evapora- 
tion in the cluster (Gieles et al. 2011). The simulations 
carried out here show that this coincides with when the 
half-mass radius reaches its maximum as well. The num- 
ber of relaxation times reached by each simulation after 
a Hubble time is given in Table 1, for those simulations 
that survive this long. Accentuated mass loss of the sim- 
ulated star clusters accelerates their dynamical evolution. 
For example, the simulated cluster evolving at Rgc = 6 
kpc undergoes 90 relaxation times in 13 Gyr compared 
to only 3 for the simulated cluster at Rgc = 100 kpc. 

5.2. Half-mass radius versus galactocentric distance at 
present time 

Figure 3 shows the half-mass radius of simulated star 
clusters orbiting at different galactocentric distances af- 
ter 13 Gyr of evolution. The values used to create Fig- 
ure 3 can be found in Table 1. After a Hubble time 
of evolution the half-mass radius defines a relationship 
with galactocentric distance which takes the mathemat- 
ical form of a hyperbolic tangent. This relation, plotted 
in Figure 3 as a solid line, is: 

r hm = rZl x t&nh(aR GC ) (5) 

where r^m is the half-mass radius and r™£ x is the 
maximum half-mass radius attained by simulated clus- 
ters evolving at 50 and 100 kpc from the galactic center 
after a Hubble time- all radii in 3D. The maximum half- 
mass radius defines a plateau in the relationship between 
half-mass radius and galactocentric distance. Current re- 
lations defining the scale sizes of clusters as a function 
of galactoce ntric distance such as the empirical power- 
law r/j = \/Rgc (van den Bergh et al. 1991) do not 
include this flattening at large galactocentric distances. 
Note that the database used by van den Bergh et al. 
(1991) only included star clusters out to 32.8 kpc from 
the galactic center. 

The parameter a is a positive coefficient that defines 
the inner slope, i.e. within the inner 20 kpc. In this case 
its numerical value is 0.06. The parameter a is a proxy 
of the tidal field of the host galaxy that is in turn due to 
galaxy mass. Note that the onset of the plateau of Figure 
3 at Rgc — 40 kpc coincides with the approximate extent 
of the disk in our models. Beyond this distance globu- 
lar clusters while remaining fully bound to the galaxy 



evolve in virtual isolation and are exempt from the trun- 
cating effects of the host galaxy tidal field. Therefore 
star clusters at 40 kpc and beyond have their sizes de- 
termined primarily by their internal dynamics (Spitzer 
& Thuan 1972) and also to some extent by their initial 
size. Hurley & Mackey (2010) showed that for clusters 
in a weak tidal field, differences in initial size can lead 
to long-term differences in half-mass radius, although the 
effect diminishes with age. Initial sizes of clusters are un- 
certain with some suggestions that the typical size-scale 
of a protocluster is ~ 1 pc (Harris et al. 2010) but the ac- 
tual size after gas removal (when the ./V-body simulations 
start) will depend on the star formation efficiency within 
the protocluster (Baumgardt & Kroupa 2007). Thus it 
should be noted that the initial size of the cluster will 
have some bearing on the location of the plateau, i.e. 

„max 
'km ' 

A set of simulations with a different initial mass was 
executed in order to establish a scaling relation that is 
independent of mass and will thus allow us to make in- 
ferences towards more massive systems. Three full sim- 
ulations, all of them with identical initial conditions but 
with different initial number of particles and thus dif- 
ferent initial masses were carried out at a galactocentric 
distance of Rgc — 8 kpc. These three simulations have 
initial masses of M = 6.3 x 10 4 M Q , M = 4.9 x 1O 4 M , 
and M = 3.2 x 10 4 M Q with initial particle numbers of 
N=100 000, 75 000, and 50 000 respectively. A unit- 
free relation that is interchangeable for these three sim- 
ulations with different initial masses is M/Mq vs t/t r h, 
where Mo is the initial cluster mass and t r h is the half- 
mass relaxation time (see also Baumgardt 2001). 

The relation above allows us to derive a predicted value 
of the half-mass radii after a Hubble time of evolution for 
simulations with twice the number of initial particles (i.e. 
N= 200 000). More massive clusters have longer relax- 
ation times t r h (Spitzer 1987). In fact, at any particular 
time, a star cluster with N= 200 000 has a t/t r h smaller 
by a factor of 3/4 than the t/t r h of a star cluster with 
half the number of particles. We predict that after a 
Hubble time of evolution, the half-mass radii of clusters 
with N— 200 000 are equivalent to the half-mass radii of 
the simulations carried out here but at an earlier stage 
of their evolution. The predicted half-mass radius values 
for models starting with A^=200 000 are plotted in Figure 
3 as blue stars. As this shows, the relation of half-mass 
radius vs. galactocentric distance postulated above (Eq. 
5) holds true for models with higher masses. 

Hwang et al. (2011) derived the two-dimensional half- 
light radii of star clusters and extended star clusters in 
the halo of the dwarf galaxy NGC 6822. The main body 
of this dwarf galaxy has a scale size of 2.3 kpc (Billett et 
al. 2002). The shape of the distribution of the effective 
radii of globular clusters from the center of NGC 6822 is 
strikingly similar to the half-mass radius vs. galactocen- 
tric distance relation shown in Figure 3. The distance at 
which the relation between half-light radius and galac- 
tocentric distance starts to plateau is about 4.5 kpc for 
NGC 6822. 

In dwarf galaxies this r™^ can be observed due to the 
smaller scales involved that can fit within the field of view 
of EST detectors. This raises the possibility of using the 
correlation between star cluster size and galactocentric 
radius to determine the structural parameters and phys- 
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ical characteristics of the different components of the host 
galaxy such as bulge mass, or disk scale-length. 

The coincidence between Fig. 3 and the results of 
Hwang et al. (2011) is encouraging, however, no infor- 
mation is available on the orbits of these clusters. A set 
of models taking into account the characteristics of NGC 
6822 will be needed before definite conclusions regarding 
the spatial distribution of the size of star clusters can be 
made. 

5.3. Core Radius 

The galactocentric distance of a star cluster also has an 
impact on its core radius, the onset of core collapse, and 
thus the ratio of core to half-mass radius. Figure 4 shows 
the evolution of the ratio of core to half-mass radius over 
15 Gyr for simulations at a galactocentric distance of 4, 
6, 10 and 100 kpc. Only selected simulations are plotted 
to preserve the clarity of the figure. The ratio r c /rh m has 
been used in various ways to characterize the dynamical 
state of globular clusters, including as a possible indica- 
tor of the presence of an intermediate-mass black hole 
(see Hurley (2007) for a discussion). 

The core radius, not directly depicted here, expands 
rapidly within the first 2 Gyr from ~ 2.6 pc to ~ 4.1 
pc. While the expansion of the core is similar at differ- 
ent galactocentric distances, the core to half-mass radius 
ratio is different for simulations evolving within and be- 
yond the inner 10 kpc of the galaxy. For simulations at 4 
and 6 kpc the ratio r c /r} lm increases significantly during 
the initial 5 Gyr of evolution. Accentuated mass loss at 
small galactocentric distances also precipitates the on- 
set of core collapse as seen in Figure 4. Core collapse is 
reached at ~ 7 Gyr for a star cluster at 4 kpc from the 
galactic center, just prior to dissolution, while core col- 
lapse is reached at 11 Gyr for a cluster at 6 kpc. Figure 
4 shows that clusters in orbits of 10 and 100 kpc do not 
reach the end of the core collapse phase in a Hubble time. 
This will continue to be true for more massive clusters 
as the relaxation time will be even larger. 

5.4. The impact of orbital ellipticity 

As a test of the impact of orbital ellipticity on the size 
of star clusters a simulation with a non-circular orbit was 
executed. This simulated star cluster has a perigalacti- 
con of 4 kpc and an apogalacticon of 8 kpc, thus inter- 
posed between the two circular orbits at 4 and 8 kpc. 
The mass-loss rate of the simulated star cluster evolv- 
ing on an elliptical orbit closely resembles the mass-loss 
rate of the cluster evolving on a circular orbit at 8 kpc. 
In fact, after a Hubble time of evolution the difference in 
mass between these two simulated clusters is only ~ 10%. 
That is, the simulation on a circular orbit at 8 kpc has 
~ 10% more mass than the simulation on an elliptical 
orbit. Similarly, the half-mass radius is ~ 12% larger for 
the circular simulation at 8 kpc after a Hubble time. 

King (1962) postulated that the tidal radius is imposed 
at perigalacticon however this has been recently debated 
(Kiipper et al. 2010b, Baumgardt et al. 2010). The ex- 
ploratory findings presented above show that the half- 
mass radius is closer to being set at apogalacticon. 

5.5. Extended outer halo star clusters 

Mackey & van den Bergh (2005) find a deficit of com- 
pact galactic young halo clusters at Rgc > 40 kpc. With 
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Fig. 3. — 3D half-mass radius vs. galactocentric distance of sim- 
ulated star clusters. Black dots give the half-mass radius of models 
at 13 Gyr. The solid line is simply a hyperbolic tangent describing 
a new relation between half-mass vs. galactocentric distance. Blue 
stars depict the predicted value of the half-mass radius after a Hub- 
ble time for simulations with twice the initial number of particles 
and twice the initial mass. 

all clusters in this region having a larger than average ef- 
fective radius, i.e. r/, > 10 pc, in agreement with the 
proposed relation of effective radius with galactocentric 
distance for star clusters 77, cx \JRgc (van den Bergh et 
al. 1991, McLaughlin 2000). 

These extended clusters are believed to be accreted 
from now-destroyed satellite dwarf galaxies with milder 
tidal fields (Mackey & van den Bergh 2005). Under the 
initial conditions of the simulations carried out in this 
study, Figure 2 shows that there is no need for merger 
events with dwarf galaxies to grow extended star clusters 
in the Milky Way at large galactocentric distances. 

The simulated star clusters evolving at 50 and 100 kpc 
from the galactic center undergo the initial expansion 
due to mass loss triggered by stellar evolution and then 
experience a small yet steady increase in half-mass radii 
(Figure 2). The tidal field is too weak at large distances 
from the galactic center to exert its truncating effects. 
Smaller initial sizes would lead to slightly smaller ob- 
served sizes. However, increased initial masses would 
more than compensate for this effect and still produce 
extended clusters at large galactocentric distances. 

5.6. A bimodal size distribution of star clusters orbiting 
dwarf galaxies 

Da Costa et al. (2009) find a bimodal distribution in 
the effective radius of globular clusters that are satellites 
of dwarf galaxies. Two peaks are seen in the size dis- 
tribution of star clusters: one at ~ 3 pc and a second 
peak at ~ 10 pc. Da Costa et al. (2009) postulate that 
the second peak at 10 pc is present when globular clus- 
ter systems evolve in a weak tidal field. They also note 
that extended star clusters in the Andromeda galaxy, the 
Large Magellanic Cloud, and the Milky Way are found 
at large galactocentric distances. 

From the .ZV-body models carried out for this work it 
can be seen that the first peak of the size distribution 
found by Da Costa et al. (2009) at r/j ~ 3 pc can be 
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Fig. 4. — Time evolution of the ratio of core over half-mass radius. 
Solid lines are models at 4 and 6 kpc while dotted lines are models 
at 10 and 100 kpc from the galactic center. Both clusters with 
orbits at 4 and 6 kpc from the galactic center reach core collapse 
before a Hubble time. 

explained by clusters that are originally massive enough 
to survive the mass loss due to tidal stripping. Their size 
however is determined by tidal interactions with the host 
galaxy. The second peak in the distribution at Th ~ 10 
pc would be created by star clusters evolving in a benign 
tidal environment, determined only by stellar evolution 
and internal dynamics. 

6. THE DISSOLUTION OF A STAR CLUSTER 

Figures 1 and 2 show that the model cluster evolving on 
a circular orbit with a radius of 4 kpc around the galac- 
tic center does not survive a Hubble time. The fierce 
mass loss of a cluster evolving at 4 kpc from the galac- 
tic center drives the cluster to dissolution after 8.4 Gyr. 
This particular simulation was run twice with a different 
seed for the initial particle distribution to ensure that 
the dissolution before a Hubble time was reproducible. 

A simulation with a more dense and compact star clus- 
ter evolving at 4 kpc from the galactic center was also 
computed. This simulation has an initial half-mass ra- 
dius of 3.1 pc and the same mass of the other simulations 
presented here, i.e. M to t ~ 6.3 x 10 4 M Q . A higher den- 
sity provides shielding from the tidal field evident from 
this denser simulation surviving up 13.3 Gyr. However 
it has less than one thousand stars remaining after 12.5 
Gyr and is only left with M tot — 764M Q . This work pro- 
vides a solid lower limit to the initial mass of star clusters 
seen today evolving at close range from the center of a 
Milky Way type galaxy. 

The dissolution of the simulated star cluster at 4 kpc 
is real and of significance to explain the galactic globular 
cluster system radial density profile. The radial density 
distribution of Galactic globular clusters is described by 
a power law with a core (Djorgovski & Meylan 1994). 
A simple power law is a good fit for the radial distribu- 
tion of globular clusters at large galactocentric distances. 
However, at Rgc < 3.5 kpc, a flattening of the radial 
distribution of clusters occurs making the overall distri- 
bution better described by a Sersic profile (Sersic, 1968) 
with a Sersic index of n=3 (Bica et al. 2006). 



TABLE 1 

Parameters of Simulated Star Clusters after a 
Hubble Time 



Rgc 




Mass 


N stars 


N binaries 


t/t rh 


(kpc) 


(PC) 


M 






4 
6 


3.5 


3.1 xlO 3 


4653 


453 


90.4 


8 


6.0 


1.3 xlO 4 


25789 


1703 


17.9 


8.5 


6.7 


1.5 xlO 4 


31602 


1916 


13.8 


10 


7.9 


1.9 xlO 4 


44393 


2331 


9.1 


20 


11.2 


3.0 xlO 4 


79591 


3575 


4.0 


50 


12.9 


3.5 xlO 4 


94666 


4069 


3.0 


100 


13.1 


3.6 xlO 4 


97506 


4161 


2.8 



Note. — Parameters of simulated star clusters after 
a Hubble time. The simulated star cluster evolving on 
a 4 kpc orbit dissolves before a Hubble time. Column 1 
gives the galactocentric distance of the circular orbit in 
which the clusters evolve; Column 2: 3D half-mass radius; 
Column 3: mass in solar masses; Column 4: number of 
stars; Column 5: number of binaries; Column 6: number 
of half-mass relaxation times that have elapsed in a Hubble 
time (13 Gyr/t r ^). 

The central flattening of the globular cluster density 
profile has been explained result of a primordial 

density distribution with a depleted core (Parmentier & 
Grebel 2005), incompleteness due to obscuration, or the 
result of enhanced destruction rates of globular clusters 
in the central regions of the galaxy. Near infrared surveys 
of the galactic core have revealed few globular clusters 
during the last two decades (Dutra & Bica 2000). In- 
completeness is thus a minor factor when accounting for 
the apparent depletion of star clusters in the inner 3.5 
kpc of the galaxy. Given that the simulation of a star 
cluster orbiting at 4 kpc from the Galaxy core shows its 
complete dissolution before it reaches a Hubble time, this 
is an argument in favor of tidal disruption being respon- 
sible for an enhanced destruction rate of star clusters in 
the inner 4 kpc of the Galaxy, particularly those at the 
lower end of the globular cluster mass function. 

7. FINAL REMARKS AND FUTURE WORK 

The parametrization of the size scale of star clusters 
presented above can be used as a primary test of the ra- 
dial distribution of extragalactic globular clusters with 
respect to the host galaxy in observational studies where 
spectroscopic information is not available. Globular clus- 
ters with r/j > 10 pc are expected to be at large gala- 
cocentric distances, extended clusters at Rgc J$ 40 kpc 
can be expected to be artifacts of projection. A well 
characterized size distribution of globular clusters across 
bulge, disk, and halo can also be used as an independent 
test of the mass of the different structural components of 
galaxies. 

The simulations with a realistic galaxy potential pre- 
sented here yield, after a Hubble time of evolution, star 
clusters with the characteristics observed today. Further 
exploring the initial values used for the set-up of the sim- 
ulations is a natural follow-up to this work. The initial 
values for the effective radius and concentration parame- 
ters can have an impact on the observed size distribution 
function (Harris et al. 2010). We will also aim to make 
a comparison with the prescriptions of star cluster evo- 
lution put forward by Alexander & Gieles (2012). 

New computing capabilities enabled by GPUs will 
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allow simulations with a larger initial number of 
star-particles. This should demonstrate that massive 
clusters can survive at small galactocentric distances, 
i.e. Rqc = 4 kpc. Such simulations will improve our 
understanding of the initial size and mass distribution 
of globular clusters and test the theories that claim that 
globular clusters are the remnant nuclei of disk galaxies 
(B6ker 2008). 
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